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Abstract 

This paper is devoted to investigation of the quantum Friedman-Robertson- Walker universe 
with matter satisfying the equation of state p = wp, where w is an almost arbitrary constant. The 
procedure starts with a reduced Lagrangian, which describes the system in a gauge fixed, so that 
the evolution parameter corresponds to the cosmological time. Then we construct the phase space, 
which is believed to correspond to the reduced phase space consisting of Dirac's observables. The 
physically relevant quantities are mapped into operators. We show that the operators have self- 
adjoint realizations and that there exist quantum states for which the evolution across singularity 
is well-defined. 
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I. INTRODUCTION 



For a gauge system like general relativity, there exist two ways of quantizing, (i) The 
Dirac method |l( starts with quantization of all the degrees of freedom of the configuration 
space of e.g. the Hilbert-Einstein action. Then the gauge freedom is removed from the 
quantum theory by solving the constraint operator equation. On the other hand, (ii) the 
reduced phase space method (see e.g. 0]) begins with canonical formulation of the classical 
theory and then the gauge freedom is removed by identification of all Dirac's observables 
through the following weak equality: {C, Q} ~ 0, where C is a constraint, Q is a Dirac's 
observable and the equality is weak in the sense that it holds on the constraint surface 
C = (it excludes the solutions of the type Q = C). Then the Dirac observables, which are 
physical degrees of freedom, are mapped into quantum operators. 

In this paper we take yet another route: we remove the gauge freedom already at the level 
of Lagrangian. Then we move to a Hamiltonian formulation and subsequently we quantize 
the canonical system using the Schrodinger representation. The analysis is restricted to the 
compact, flat FRW universe. 

In section [Til starting with a given Lagrangian, we show that it describes the system 
under consideration and arrive at canonical formulation in convenient basic variables. In 
section [LTTl we map the phase space functions into operators and study their self-adjointness. 
Next, in section [IV] we compute the evolution of mean values of the relevant operators. We 
conclude in section [V] 

In appendix [X] we make conjecture that the reduced phase space method is closely related 
to the reduced Lagrangian method. 

In appendix [B] we report on an unsuccessful attempt at resolving the singularity. It is 
discussed and compared with approach used throughout this paper in concluding section [V] 



II. CLASSICAL THEORY 

The metric of the flat FRW universe in the comoving coordinates and with respect to the 
cosmological time t reads: 



where a(t) is a scale factor and dx is a distance measure in a space-like leaf E. Instead of 
the scale factor a, we are going to consider the dynamics of the physical length between 
two unspecified points in £ and we denote it by l{t) := h-^pj, where 1$ is the length at the 
moment t = to- Therefore, for each value of t, the length l(t) is a physical quantity, i.e. the 
Dirac observable. Now let us see how one can introduce the dynamics of /. 

A. Lagrangian formulation 

Let us examine the following action integral: 



ds 2 = -dt 2 + a(t) 2 dx 2 , 



(1) 




(2) 
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where a := • vg^ and A is a constant of the length dimension and will be specified later. 
The parameter w is a real constant. Variation with respect to I gives the equation of motion: 

- (3w + l)al 3w l 2 - 2al 3w+1 l = = (/^) , u , (3) 



which has the solutions: 



3w + 3 I 



l(t) = lo[^ T ^f o (t-t ) + lj . (4) 
satisfying the initial conditions: 

/(to) = lo, i(t ) = k. (5) 
Let us turn to the Friedman equations: 

H = — p, (6) 

H 2 + H=-^(p + 3p), (7) 

where H is Hubble's parameter, p is the energy density and p is the pressure of the matter 
in the universe. Now, making use of H := I, we find out that: 

# 2 = ^o W+1 )j^r 3 , H* + H=- 1 -(l + S^/r 1 )^' (8) 

and we easily conclude that our system given by the action integral ([2]) is the flat FRW 
universe with ^f-p = (/'o/o W+1 )i3^+3 and p = wp. 

Hence, we have succeeded in the gauge-fixed Lagrangian formulation of the dynamics 
of the flat FRW universe filled with matter satisfying the equation of state p = wp. The 
connection of this reduced Lagrangian with the reduced phase space procedure, which begins 
with the introduction of the kinematical phase space (consisting of both gauge and physical 
degrees of freedom) , is commented in appendix [A] 

The energy of matter contained in a fiducial cubic cell of the edge length changing with 
time as l(t), reads: 

3/ 3 // \V/ X ~ 3 ™ 



and, for w ^ 0, is clearly not conserved. For instance, the matter may consist of particles, 
which move with respect to the comoving reference frame, then the particles loose/gain 
their kinetic energy as the spacetime expands/shrinks and the energy of matter changes 
with the evolution of the universe. However, one may introduce a conserved quantity, the 
total energy: 

E *ot = ■ / 3W+3 , (10) 

which is equal to the energy of matter at the moment when the value of the length of the 
fiducial cubic cell is equal to A. If the moment happens to be t = t , then E tot = E matt (to) = 

||(|) 2 ,andA = / . 



1 Note that these equations refer to the specific choice of evolution parameter, which is the cosmological 
time, and hence there are no gauge degrees of freedom in this formulation. 
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B. Canonical formulation 



We perform the Legendre transformation: 

I i-> I e R+, i^ Pl :=—r = 2al 3w+ H E R, (11) 

dZ 

so that 

H = p / /-L = i -^ TT and {l, Pl } = l (12) 



Hamilton's equations: 



Pi ■ r ran to , i\ 



/ = {/ ' H} = 2^i' ^ = te^} = ( 3 ^ + 1 )4^2- (13) 

can be integrated and we obtain: 

= CI 2+ , or = 0, I = C (stationary universe), (14) 

where C is any real constant. The solutions split into three cases: w < —1/3, w = —1/3 
and w > —1/3. They respectively correspond to the momentum of space blowing up, being 
constant and vanishing as the universe approaches the singularity 1 = 0. 

We are going to study physically interesting quantities like Hubble's parameter or energy 
density, which take the following form, respectively: 

The Hamiltonian EI is a conserved quantity like the total energy of the system defined in 
eq. (IT0|) . The comparison of these two leads to the condition: 

1 Pi M=F J_J P± L n 6] 

aUto+i M ^ a2 87rG4/ 3„ + i A 3«,' ^ i0 ^ 

and hence 

01 = 8^' (17) 
which is in accordance with our prior definition and explains the physical meaning of A. 



C. Canonical transformation 

In what follows we will look for a more convenient conjugate pair. We assume that: 

1. the canonical transformation is valid globally 

2. the 'stationary universes' in ( |T4l) are included in the system and thus well-mapped by 
the transformation 

3. canonical transformation does not extend the phase space, in particular it preserves 
the topology of the phase space under consideration 
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Now, let us break the above rules in the following example: 

{l, Pl } — > {1 := H,p ; := 3^^ } , (18) 

The map is singular for pi = and thus the 'stationary universes' being the points at the 
axis pi — 0, are excluded now. Instead, the closure of the range of the map covers a new 
region of the space of parameters, p~\ = 0. As a consequence, the trajectories in the phase 
space are complete and the singularities are absent! But this is not a good way to resolve 
the singularity, since, apart from the lack of quantum physics, the procedure has a drawback 
in that it is non-unique and here are further examples: 

{l, Pl } — > {! :=M,p l : = J^-^ + sgn(p ; )/(H)} , w ^ -1, /arbitrary. (19) 

Therefore, making use of a 'regular' map, we introduce the following new canonical pair: 

{; - pi} - 0" : = ^r/" +3) ' 2 ^ : = 2i^ ] • w * ~ L (20) 



Now, the Hamiltonian reads: 
the Hubble parameter: 
the energy density: 
and the length: 



m = -pi (2i) 



a 



H= l --^—K (22) 

a 3w + 3 I v ' 

o 1 2 g ^ 

P _ a* vrG(3 W + 1) [2 ' ^ 

z= ( 3^ + 3 r) 2/(3. + 3) ^ 



III. SCHRODINGER QUANTIZATION 

In the Schrodinger representation: 

I ^ 1 ■= x, pi ^ p l ■= -ih-^, (25) 

ax 

and we get the following formally self-adjoint operators that are relevant in the cosmological 
context and will be studied below: 

ti 2 d 2 

H^H=---f- (26) 
a dx A 

H H = — - ~~~^(~~ t ~~j~ x - ™~ V , (27) 
a 3u> + 3 ax ax 

P ^ ^ = ^ 2 7 rG(3l + l) ( -^" 1 - ^"h^ (28) 

/^/=(^±^x) 2/(3w+3) . (29) 

The range of x is (0, oo). Thus, the problem of evolution across singularity is formulated in 
terms of a freely moving quantum particle on a half-line. 
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A. Hamiltonian and the evolution 



The theory of the operator — 4^ on half-line may be found in [3]. In short, there are 

infinitely many ways, in which one may define the operator EI to be essentially self-adjoint on 
L 2 (R + , dx) and we will denote these options by H M for fi e R U {oo}d. In order to calculate 
the evolution of a given state one may apply the appropriate Fourier Transform, denoted by 
FT^, and multiply the resultant wave-function by e~ lk t , where k is the momentum. 

The important point to be made here is that we do introduce an evolution operator. In the 
quantization scheme proposed in [2], one only quantizes the Dirac observables enumerated 
by a classical evolution parameter - one does not introduce a generator of change in this 
parameter, called the true Hamiltonian, in order to quantize it. However, in our model 
of universe this approach would not work - it would not give any sensible solution to the 
singularity problem. The crucial difference is that in the model presented in [2j the classical 
dynamics was already non-singular, whereas in the model considered here the dynamics ends 
at the point I = 0. 

The quantization of the evolution may give some hopes to resolve the singularity. The 
following heuristic argument shows that the closer to the singularity the more the quantum 
universe departures from the classical one: 

{H,U} ~ -2H 2 , (30) 

\[H, H] ~ —2H 2 + -x~ 4 . (31) 
% 2 

We observe that at the classical level (I30p the value of the Hubble parameter H goes to 
minus infinity as the universe approaches the singularity, whereas at the quantum level (IBTj) 
the term x _4 opposes the unbounded decrease of the operator H. From the form of H given 
in (127)1 . one sees that the 'closer' the wave-function to x = 0, the bigger the value of H 2 . 
But at the same time, the closer the wave-function to x = 0, the bigger the value of the 
opposing term x~ 4 . 



B. The energy density and Hubble operators 

We will study the formally self-adjoint Hubble operator: 

3w + 3 a ^ dx dx^ ^ ^ 

Let us consider the following inverse mapping L 2 (R + ,dx) h-> L 2 (M + , dy): 

$(x) ^ <j>{y) := y>0. (33) 

Let us see that this mapping is isometric, and hence unitary: 

^(xMx) dx = r Ii{2y/V)1>(2>/V) ^ = f° } dy. (34) 

o Jo Jo y y 



Usually, the choices fj, = and fi = oo correspond to the Neumann and Dirichlet conditions, respectively. 
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Application of the mapping ( 1531) to the operator H gives: 



^ -y~ 1/4 (-x- 2 + 2x-^)y^ = (35) 
Sw + 3 a ax Sw + Saay 

which is a simple derivative operator. Analogically, under the same mapping the energy 
density operator p is transformed into Laplacian ~ — J^-. 

Again, there are infinitely many inequivalent ways, in which one may define a self-adjoint 
p, and we will denote these options by p^, p E f U {oo}. However, the Hubble operator, 
being a momentum operator, cannot be self-adjoint on the half-line. Therefore, we will 
redefine H := yfp. The spectrum of these operators reads sp(p^) = sp(Jf M ) = M + for any p. 



C. The volume operator 

The volume operator v := I s = ( 3w + 3 x ) 2 / ( - ul+1 - ) is a self-adjoint operator on L 2 (R +J dx) 
and its spectrum for w > — 1 reads sp(v) = R + . 



IV. EVOLUTION OF PHYSICAL QUANTITIES 

Let us restrict to the study of the evolution that preserves the Dirichlet condition (i.e. 
generated by for p = oo). We pick the following family of wave-functions: 

4> a (k) := Nke~ ak '\ k > 0, (36) 

where N = 2 ^ 2o jy 4 , [A;] = [x]^ 1 = m~^2^, [a] = m 3w+3 and evolve them in time: 

exp (-i*H)0 o (Jfe) = exp (-i-tk 2 )j> a (k) = Nke- {a+i ^ t)k2 = $ a (k, t). (37) 



Applying the Fourier sine transform FT 00 renders: 



Nx 



4> a (x,t) := / Bm(kx)<f> a (k,t)dk = — —w^Jl e . (38) 

Jo (A a + i>-t)) 1 



A. The invariant w+ \fE^ 

Let us calculate the expectation value of the operator corresponding to the classical 
invariant w+ -\/E tot : 

( w+ <fEZt ) ■= (k{k)\ "VI Mk)) = w+ \f^\N\* f e-^k^ dk 

V a J 

w+1 [h 2 ( 2 /5 + 3uA\ w+1 1 



2a\ v ^ T \2 + 2w) ) 'a' ^ 
It will turn out to be useful for linking the classical trajectories with quantum states </>„. 
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B. Volume 



The expectation value for the volume operator for the state <f) a (x,t) equals: 

/ u ri N2 (^ + 3\ 2/(w+1) r , 

(v)(t) = to : — / x ™+ 1 e ^ ax — 

N 2 {3w + 3\ 2/(w+1) i A i / a \-i-iKY3 

+ - — 



/3u. + 3V^".^/ a V r( , + _, (40) 



T )a 2 + —t 2 )^. 



a J \2 + 2wr a 2 



One observes that at t = the volume acquires its minimal value 



w-2 

i 2™+i , 2 Tn /5 + 3to , 

u^:=a^— ^(3w + 3)^r — — . (41) 
WTT \2 + 2w' 



Let us study the asymptotic behavior in (j4"0l) for i — >■ ±oo: 

WW-^(3u. + 3)Arf|±^V5±l , 1C>^. (42) 



On the other hand, the classical solution in (J4J) tells that the volume of the fiducial cubic 
cell with the edge I evolves as: 

«t) = *(^{t- toft + (43) 



2 v u/ / 



which for t — t ^> 3^3^ reads: 



2 , 2 



U)+l 



^)-^(^)"\3^ + 3)^(t 2 )^= -y^(^) OT. (44) 

The comparison of (jUJ) with (I42p gives the following relation: 

h 2 ( 2 /5 + 3™\V +1 1 



2a\^ V \2 + 2w)) E tot ' ^ 



which is in accordance with ( 1551) . if we substitute E tot for ( "yE^)™" 1-1 . Plugging this 
relation into the formula (141 p we find out that at the bounce the critical volume of the cubic 
cell reads: 

v crit = V ( 5 -±^) (*0±1) ^ -(/^ . (46) 

First, notice that v cr a depends only on the classical invariant of motion E to t, so it does not 
depend on the particular moment of evolution of the classical universe, at which we calculate 
v cr it- Next, notice that it is proportional to — , hence the more massive the universe the 

smaller its critical sizeu- The formula (T46|) knows nothing about Planck scale and combined 
with the knowledge of the early Universe evolution, can lead to the bounds on the mass of 
the whole Universe. 



3 The more massive universe means the universe of bigger size with the same energy density of matter or 
the universe of the same size with higher energy density of matter. 
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FIG. 1. The comparison of classical and quantum behavior of the volume for w = 1/3, a = H = 
a = l. 

C. Energy density and Hubble parameter 

First we make use of the unitary mapping L 2 {R + ,dx) L 2 (R + ,dy) and map the state 
4> a {x,t) accordingly to (133|) . We find: 

Nv 1/4 y 

4> a (x,t) ^ Mv,t) ■= y-^V-CVy,*) = (2(a + ^ ))3/2 e ~ w - ( 47 ) 

We assume that the action of H and p on ip(y) G L 2 (W + , dy) preserves the Dirichlet condition 
(i.e. we set p = poo). We perform the Fourier sine transform FT 00 on ip a (y,t): 

~ , f°° , s ~ sin (| arctan((a + it)k)) , 

Mk, t) = / sin(kx)A(x, t)dk = N U 1 " ' >\ (48) 



where N = ^^3/2 an d hence we have 



Ma = \N\* 



sin(| arctan(^)) 



(C 2 + k 2 )i 



(49) 



1 A 

where C = , _^ if , . Now, in the momentum representation, the operators H and p are 
proportional to k and k 2 , respectively. Since the function ( 1491) for large /c behaves like 



$J a ~ A;" 5 / 2 , (50) 

only the operators O such that O < Ck a for large k, a constant C and a < |, are well 
defined on the states that we consider in the sense that: 

WalOtPa) < 00. (51) 

However, we can still consider the expectation values 

$a\d$ a ) ■■= {OHalOHa), (52) 



9 



H 



3.0 






£.5 






1 2.0 






/ 15 






/ 1.0 




0.5 









' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — I'll t 

-10 -5 5 10 



FIG. 2. The comparison of classical and quantum behavior of the Hubble parameter for w = 1/3, 
a = H = a = 1. Just after the Big Bounce the quantum universe is shown to expand faster than 
the classical one, (H) > H. Thus, the quantum universe seems to undergo a kind of inflationary 
phase. 



for operators O < Ck . The reason is that the integral (ip a \k tp a ) is convergent and since 
the states ip a may be approximated by states from the domain of the operator k 2a up to 
arbitrary accuracy, the integral ( 1521) keeps the physical meaning of the expectation value of 
this operator. 

In particular, we may obtain the integral formula for the expectation value of the Hubble 
parameter (but not the energy density): 



(i>a\k$ a ) = \N\ 



sin(| arctan(^)) 



(C 2 + k 2 )\ 



k dk < oo. 



(53) 



We have already learnt from the study of the volume that the bounce occurs for the states 
ip a at t = 0, hence we may calculate the expectation value of the operator H = 3^3^ at 
the bounce: 



H, 



crit 



h 

a 



4\N\' 



-<22 



3w + 3 

After plugging the value of a from 



sin(| arctan(fc)) 



k dk = — 

a 



8T 



4' 



/ 27T 



l+k 2 )l 

51) into fl54l). we have: 



(3w + 3) 



5 
3 



1 

a 



(54) 



H, 



ier 2 (f) 



/ 27T 



crit 



(3w + 3) 



5 
3 



E, 



tot 



( 2 r ( 5+3w \\ 
Vv^F \2+2wJ) 



w+1 



h 



(55) 



One immediately sees that the value of Hubble parameter at which the universe 'reverses' is 
proportional to the total energy of this universe. Therefore, the more massive the universe, 
the higher the value of H CT u. Moreover, one may safely assume that the energy density at 
the bounce p crit is proportional to the square of the total energy 



E 2 

11 'tot- 



V. DISCUSSION 



It is a very interesting conjecture to think that our phase space with the physical quan- 
tities like the Hubble parameter and the Hamiltonian is in one to one correspondence with 
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the reduced phase space arising from the identification of Dirac's observables in the full 
canonical gauge-symmetry formulation (for a short comment see appendix [A]). If so, then 
our Hamiltonian corresponds to a true Hamiltonian in the reduced phase space formulation. 
This conjecture is investigated in [ij and a special attention is given to the freedom in fixing 
an evolution parameter. 

Our analysis shows a way to resolve the singularities of classical theory. It started with 
a Lagrangian that was free from gauge symmetries. The evolution parameter t was treated 
as a real, absolute time in which the system evolves. This cosmological clock does not exist 
outside the universe and should be connected with some measurements or change that we 
observe on cosmological scales. Nevertheless, no matter which quantization method one 
uses for a gauge system with a Hamiltonian itself being a constraint, one needs an evolution 
parameter. We argue that the choice of the evolution parameter in such a way that the 
singularity occurs at its finite value is crucial. Thus a popular free scalar field is not a good 
choice and we encourage the reader to go through appendix [B] to see this. 

In our approach we managed to express the dynamics of the FRW universe in terms of 
a freely-moving non-relativistic particle on half-line. The classical dynamics ends when the 
particle hits the boundary of half-line and the evolution parameter (i.e. the cosmological 
time) cannot be extended beyond this moment. However, we can promote the Hamiltonian 
to a self-adjoint operator in quantum theory. Due to the Stone theorem, the self-adjoint 
Hamiltonian can be exponentiated to a unitary operator, which gives the evolution of the 
system for the unbounded values of the evolution parameter. Thus, the half-line of the 
cosmological time is extended to the whole real domain and the singularity is resolved. 

We have studied the evolution of quantum states a , where a G M+. How much of the 
Hilbert space is covered by Span{0 a }? We study this question and others concerning the 
evolution of unbounded operators p and H in p . 

We observe in the evolution of the Hubble parameter and the volume of the universe 
no fixed scale, and in particular Planck scale is irrelevant. The critical values v cr u and 
H crit depend on how massive the universe is in the classical phase of its evolution, which 
corresponds to the value of a. However, we can calculate the value of (p){v 1+w ) at the 
bounce: 

h 2 /2T( 5+3w ) \ w+1 /k\ 

*Hl it v%? = • 5 > + 1) (-^F J T 2 (j J , (56) 

which is clearly independent of the parameter a, i.e. of a particular quantum state. Although 
this quantity is energy dimensionally it is not the same with the energy of the system. Notice 
that the mean value of energy of the system obtained in (|39|) depends on a state via the value 
of a. Thus, ( 156]) may signal the existence of some fundamental energy scale not connected 
with any particular universe. The existence of such scale may enable to fix the value a. 

This brings us to the problem of determining the value of a. First, notice that the choice 
of a is connected with the choice of time and to this extent it is arbitrary. One may suspect 
that the study of the conjectured correspondence (appendix [SJ may shed light on this issues 
and thus we postpone the investigation to 

Another indeterminacy of the system is connected with the fact the addition of a total 
derivative to the Lagrangian in (j2J) leads to the same dynamics. More importantly, there is 
an infinite number of ways in which one may define the self-adjoint Hamiltonian. We have 
restricted to a single choice, preserving the Dirichlet condition. But other choices seem to 
be equally good and are studied in jBf. 
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Appendix A: Relation between reduced phase space equipped with M trU e and re- 
duced Lagrangian I^true 

Suppose one has constructed the reduced phase space equipped with the induced Poisson 
bracket and a true Hamiltonian M true . In principle, one may now construct an action integral, 
that is an integral of the Lagrangian L 4rne defined as a function of the Dirac observables: 

Ltrue = P' J Oj> — EIt rne , (Al) 

where P J ' denote the canonical conjugates to Oj> (notice that the number of the Dirac 
observables is always even) and Oji = {Oj>, M. trU e\ind- 

Therefore, one may begin with a reduced form of Lagrangian L true , which encodes the 
dynamics of the system in a gauge fixed in correspondence with M true . This formulation 
should also enjoy a freedom in choice of evolution parameter, hence there should be a class 
of h true and a corresponding class of M. true . 



Appendix B: Massless scalar field case 



The metric of FRW universe reads: 

ds 2 = -N 2 (r)dr 2 + a(r) 2 dx 2 , (Bl) 

where r is an evolution parameter. The choice of the shift function N corresponds to the 
choice of the evolution parameter. The rest of the notation as well as the meaning of the 
variable I used below were explained in section HT1 



1. Lagrangian formulation 

For the flat FRW universe with massless scalar field the action reads (see e.g. p): 
which leads to the Euler-Lagrange equations: 

72 Q 72 12 



5S 3 d ( 2ll\ 3 I 2 3/ 2 

-TT = 0^— — - — --—- + -- 



51 16ttG dt\N J 16nGN 2 N 

5N NA 16ttG + 2 9 
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to which the solutions in the cosmological time (N = 1) read: 

1 



(8vrG^i)(t + to) 



\n(t + 1 ) + 



3Q-2 



(B3) 



2. Canonical formulation 

The Legendre transformation followed by Dirac's analysis gives: 

,2 - 



\ 3 I 2 I 3 J 



(B4) 



— . The Hamiltonian 



where the phase space is 4D: (l,pi, where pi = g%^, V<t> z 

EI is a constraint and N is a non-zero coefficient corresponding to the choice of gauge. 
Let us look for the Dirac observables via fixing the gauge N = I 3 and solving: 



{O i; H} = 



8ttG 7 / d 



4>+ 



'd0 



O,; = 0. 



(B5) 



We introduce the new 'geometrical' variables X± = ha. pi ± In/ so that the equation (1B5I) 
reads now: 

ri R^n 



16nG v d d 



dx _ +P* d(j) 



Oi = 0, 



and the solutions are easily found to be: 

o 1 = x + , 2 = i 



2 



8ttG e*+ 



(B6) 



(B7) 



which are not independent on the constraint surface due to the identity HI = ^pe 2C>1 — \0\ 
and for which the Poisson bracket is now modified as follows dl A dpi = \e x+ dX + A g?X_. So 
the complete set of independent elementary Dirac observables consists of two elements, say 
0\ and 3 , satisfying the algebra: 

{0 1 ,0 3 } = e'°\ (B8) 

We redefine f2i = e Gl and VL 2 = O3 — | In so that 

{J2 1 ,fi 2 } = l, (B9) 

and we restrict the analysis to the 'expanding universe' solutions Qi = Ipi > and on the 
constraint surface we have O3 = |X_ + From (1B3|) we have that: 



L = ^87rGfiie- 3n 2e A 



3 



(BIO) 



This quantity is now the Dirac observable for a fixed value of <j), since it does not depend 
on the value of t that one ascribes to the moment of measurement. The evolution of length 
is defined as the flow through the parameter <ft, which parameterizes the family of Dirac's 
observables L. 
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3. Classical evolution 



Let us identify the generator of the change in <p: 



= ^ * + if = - H ^ + «s ln - (B11) 

where / is arbitrary. How to get rid of this indeterminacy? Note that to parameterize the 
reduced phase space one needs two independent Dirac's observables. The first choice was 1$ 
and what is to be the second? Let us obtain from (1B3j) the Hubble parameter: 

H(j> = ie^e" 3 ^^, (B12) 
3 

so that 



{H^M = dtHt * ^ = yj ^ O = yj 5^fl! + ^(na), (B13) 

where g is arbitrary. Now, comparing ( IBllj) and (1B13j) . removes the indeterminacy and the 
true Hamiltonian is := y/^pOi. Thus, instead of considering the whole family of 1$ and 

Hfj,, which of course produce redundancy in parameterizing the reduced phase space, one 
fixes the value of (f> and includes the as a dynamics generator. One only needs to include 
that: 



{H^} = ~-1^, {l+,M+} = J^l+, {H+,W+} = -3y]^H+, (B14) 



where = J ^^H^. Now, let us identify p^, the conjugate of If. 

{l+,pd = l^V<t> = (SnGy^VW*^^ + /(n ie -»») = -^5#* + f(k), (B15) 
where / is arbitrary but we put / = 0. 



4. Quantum evolution 

Now, let us collect the information on the system in the form convenient for quantization: 

{Z,p4 = l, H=W— Zp^, /£ ,= ____J j w = Z d , (B16) 

where Z > 0, we have dropped the index <p for convenience and thus Z, HI, H and p are the 
length, the Hamiltonian, the Hubble parameter and the energy density, respectively. 
We do the following canonical transformation: 

l^x = h\ p^ p = it, (B17) 

where the new canonical pair {x,p} = 1 is also defined on M. + x R and 



H = V24nGxp, H = —r-p, p = —z~P, v = 3x. (B18) 

o o 
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Now, the Schrddinger quantization, p i— > p =: — i-4- and x H- x = x, renders the following 
formally self-adjoint operators: 

d d \ 8nG ( d \ 3 



H = V6^G[ - ix— - i—x , H = - i— ), p = — -H 2 , v = 3x. (B19) 

\ ax ax J 3 \ ax J 8nG 

Let us consider the following unitary mapping: 

R + 3 x, ^(x) i — y $(y) = V2e y V(e 2y ), y e R, (B20) 
and we find that the Hamiltonian operator in this new Hilbert space reads: 

U=V6^G(-i-^-). (B21) 
dy' 

Now we should ensure that the Hamiltonian is a positive operator, but we will ignore it for it 
will not alter the fate of the singularity. Now, we can obtain the evolution of wave-functions: 

e *H*$( y ) = $( y + <pV0^G) = V2e y e (f ' V ^y{e 2y e 2(l ' V ^), (B22) 

hence the inverse mapping back to the start Hilbert space gives: 

^y^V^G^^y^V^G^ ^ ^V^G^^V^G^ ( B23 ) 

Now, we are ready to study the evolution of the expectation value of p: 

(e iH ^(x)|p e iH ^(x)> = ^ e 2v^0 f _ ^*( x ) dx, (B24) 

o j ace 

and the expectation value of v: 

(e m ^(x)\v e ffl ^(x)) = 3e- V2 ^^ > J dx. (B25) 

The classical behavior is reproduced in the quantum theory, which includes the unbounded 
growth of energy density and the vanishing size of the universe. 

One may think that this shows that the singularity cannot be resolved in the above system. 
In fact, however, the singularity problem was not addressed at all due to the peculiar choice 
of the evolution parameter 0. In order to study the structure of the singularity one has to 
ensure that the classical singularity occurs at finite value of an evolution parameter. 
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